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We point out the opportunity of the tau neutrino astronomy for the neutrino energy E ranging be- 
tween 10 GeV and 10 3 GeV. In this energy range, the intrinsic tau neutrino production is suppressed 
relative to the intrinsic muon neutrino production. Any sizable tau neutrino flux may thus arise 
because of the — > v T neutrino oscillations only. It is demonstrated that, in the presence of the 
neutrino oscillations, consideration of the neutrino flavor dependence in the background atmospheric 
neutrino flux leads to the drastically different prospects between the observation of the astrophysical 
muon neutrinos and that of the astrophysical tau neutrinos. Taking the galactic-plane neutrino flux 
as the targeted astrophysical source, we have found that the galactic-plane tau neutrino flux domi- 
nates over the atmospheric tau neutrino flux for E > 10 GeV. Hence, the galactic-plane can at least 
in principle be seen through the tau neutrinos with energies just greater than 10 GeV. In a sharp 
contrast, the galactic-plane muon neutrino flux is overwhelmed by its atmospheric background until 
E > 10 6 GeV. 

PACS numbers: 98.38.-j, 13.85.Tp, 14.60.Pq 



I. INTRODUCTION 

The v p — * v T neutrino oscillations established by the 
high-statistics Super-Kamiokande (SK) detector, ensure 
that a non-negligible v T flux reaches the Earth. A recent 
SK analysis of the atmospheric neutrino data implies the 
following range of the neutrino mixing parameters |lj 

Sm 2 = (1.9 - 3.0) • 1CT 3 eV 2 , sin 2 20 > 0.9. (1) 

This is a 90% C.L. range of the neutrino mixing parame- 
ters with the best fit values given by Sm 2 = 2.4- 10 -3 eV 2 
and sin 2 20 —1 respectively. 

The tau neutrinos resulting from the above — > v r 
oscillations are presently identified on the statistical basis 
. On the other hand, the total number of observed non- 
tau neutrinos from various detectors are already greater 
than ~ 10 4 with energies ranging from ~ 10 _1 GeV to ~ 
10 3 GeV 0] . It is essential to develop efficient techniques 
to identify the tau neutrinos jj]. 

There are at least two important reasons for observ- 
ing the v T . First, seeing the v r confirms the v p — * v r 
oscillation interpretation for the atmospheric neutrino 
data. Second, since the atmospheric v T flux is gener- 
ally suppressed as compared to the atmospheric v p flux, 
the prospective observation of the astrophysical v T suffers 
much less background than in the case. 

In this paper, we address the second point with the 
galactic-plane tau neutrinos as our illustrating astrophys- 
ical source. We point out that contrary to the general 
expectations, the atmospheric neutrino flux does not 
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dominate over the astrophysical neutrino flux for the neu- 
trino energy E < 10 3 GeV, once the flavor composition 
of the neutrino flux is taken into account. The idea for 
such an investigation has appeared earlier in Ref. 

In the context of two neutrino flavors, and u T , the 
total tau neutrino flux arriving at the detector on Earth, 
after traversing a distance L, is 

<j>™(E) = P{E) ■ <f> Uii (E) + (1 - P(E)) ■ <f>^(E), (2) 

where P{E) = P{v^ — > u T ) = sin 2 20 ■ sin 2 (L/L osc ) with 
the neutrino oscillation length given by L sc = ^E/8m 2 . 

In order to compute the total v T flux from a given 
astrophysical site, we need to first compute the intrinsic 
as well as the intrinsic v T flux from the same site. 

II. GALACTIC-PLANE NEUTRINO FLUX 

One calculates the intrinsic galactic-plane v p and v r 
fluxes by considering the collisions of incident cosmic- 
ray protons with the interstellar medium. The fluxes are 
given by 

ME) = Rn P J dE P MEp) app ^ +Y , (3) 

where E p is the energy of the incident cosmic-ray proton, 
dcTpp—tv+Y /d_E is the v energy spectrum in the pp colli- 
sions. R is the typical distance in the galaxy along the 
galactic-plane, which we take as 10 kpc (1 pc ~ 3 • 10 13 
km). The density of the interstellar medium n p along the 
galactic plane is taken to be ~ 1 proton per cm 3 . The 
primary cosmic-ray proton flux, cf) p (E p ) = dN p /dE p , is 
given by 

/ , n-2.74 

MEp) = 1-49- (E P + 2.15 ■ exp(-0.21 v^)j , (4) 
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Intrinsic Muon And Tau Neutrino Fluxes 
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FIG. 1: The intrinsic galactic-plane and the intrinsic down- 
ward going atmospheric muon and tau neutrino fluxes as a 
function of the neutrino energy in GeV. Details are provided 
in the text. 



in units of cm _2 s _1 sr _1 GeV _1 . The above flux is un- 
der the assumption that the cosmic-ray flux spectrum in 
the galaxy is a constant and equal to its locally observed 
value. The galactic-plane neutrino flux (abbreviated here 
as Gal) is sometimes also referred to as the galactic cen- 
ter region neutrino flux (abbreviated as G or GC), the 
galactic disk neutrino flux or the Milky Way neutrino 
flux. We shall estimate here the neutrino flux coming 
from the galactic-plane direction only as transverse to it, 
the n p decreases essentially exponentially |j| , and so does 
the <j> v (E) according to Eq. @. 

The neutrino production process p + p^is + Y is me- 
diated by the production and the decays of the tt, the K, 
and the charmed hadrons. The galactic muon neutrinos 
mainly come from the two-body n decays and the sub- 
sequent three-body muon decays. While the decay rates 
and the decay distributions of the tt and the [i are well 
understood, the differential cross section for the process 
p + p ^ it + Y is model dependent. We adopt the pa- 
rameterization for this cross section from Ref. , which 
is obtained by using the accelerator data in the sub-TeV 
energy range 0. We remark that our galactic-plane 
flux compares well with a previous calculation using 
the PYTHIA 

The galactic-plane v T flux arises from the production 
and the decays of the D s mesons. It has been found to 
be rather suppressed compared to the corresponding 
flux Fig. n shows the intrinsic galactic-plane 

and v T fluxes obtained by using the above description 
(along with the corresponding intrinsic downward going 



atmospheric neutrino fluxes). We shall use these in our 
subsequent estimates. 

The total galactic-plane tau neutrino flux, 0^ ot (_E), 
is therefore dominated by the — > v T oscillations in- 
dicated by the term P{E) ■ 4> Vfl [E) in Eq. ©. With 
the best-fit values for the neutrino mixing parameters, 
we have ^°^{E) ~ <j> v {E)/2, neglecting the contribu- 
tion of <j> VT {E), since L osc -C L, where L ~ 5 kpc. 
The total galactic-plane tau neutrino flux, </>*° t (£') = 
dN/d(\og 10 E), can be parameterized for 1 GeV < E < 
10 3 GeV, as 

&w =A ((m) a > (5) 

where A = 2 ■ 10~ 5 is in units of cm _2 s _1 sr _1 with a = 
-1.64. 



III. ATMOSPHERIC NEUTRINO FLUX 

The calculation of the total atmospheric tau neutrino 
flux is more involved. We follow the approach in Ref. 
@ for computing the intrinsic atmospheric flux which 
can oscillate into v T . 

For the 7r-decay contribution, the flux formula reads: 

aMx - J E dEN J E ^dMi-i.) 

f x dX' 

Jo *N 

x—F Nn (E^,E N ) 

x exp 4>n{E n ), (6) 

where E is the neutrino energy and £ is the zenith an- 
gle in the direction of the incident cosmic-ray nucle- 
ons. The j K = m^/m* and is the pion decay length 
in units of g/cm 2 . The A at is the nucleon interaction 
length while A^r is the corresponding nucleon attenua- 
tion length. 0Ar(i?Ar) is the primary cosmic-ray spec- 
trum. We only consider the proton component of 4>n, 
which is given by Eq. We take A at = 86 g/cm 2 and 
A at = 120 g/cm 2 0. The function P 7r (E 7I ,X, X') is the 
probability that a charged pion produced at the slant 
depth X' (g/cm 2 ) survives to the slant depth X (> X'). 
The Fjs/ 7T (Et T , En) is the normalized inclusive cross sec- 
tion for N + air — » 7r ± + Y, and is given in the Ref. [g. 

The kaon contribution to the atmospheric flux has 
the same form as Eq. (0 with an inclusion of the branch- 
ing ratio B{K — ► [w) = 0.635 and appropriate replace- 
ments in the kinematics factors as well as in the normal- 
ized inclusive cross section. We remark that the current 
approach neglects the 3-body muon-decay contribution 
to the flux. This is a good approximation for E > 10 
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The Fractional Contributions 
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FIG. 2: The fraction of contributions by the 7r, the K, and 
the charm decays to the overall downward going atmospheric 
z/ M flux as a function of the neutrino energy in GeV. 



GeV [Sj. The relevance of the muon-decay contribution 
for 1 GeV < E < 10 GeV will be commented later. 

We stress that the ir and the K decays are not the only 
sources for the atmospheric flux. For E > 10 6 GeV, 
the charm-decay contribution becomes more important 
than those of the ir and the K decays. We have used 
the results from the perturbative QCD to estimate this 
contribution [T^ . The muon neutrino flux due to charm 
contribution can be written as 



d 2 A^ 



(E,X) 



ZphZhv„ exp(-X/A p )(/) p (E) 



dEdX 



^ 1 

h 



Z P p(E) 



(7) 

where h stands for the D^, the D°, the D s and the A c 
hadrons. The Z moments on the RHS of the equation 
are defined by 



Zij(Ej) 



dE„ 



4>i{Ei) Xi{Ej) dniA^jviEiiEj) 



dE» 



with dmA^jY(Ei,Ej) = da lA ^]Y{Ei,E j )/aiA{Ei). 



(8) 
In 



the decay process, the scattering length is replaced by 
the decay length di , whereas the driiA^jY (Ei , Ej ) is re- 
placed by the dT^jy (Ei, Ej)/Ti(Ei). We note that this 
part of the atmospheric flux is isotropic, unlike the 
contributions from the 7r and the K decays. This differ- 
ence is attributable to the lifetime difference between the 
charm and the -k(K) mesons. The various Z moments are 
calculated using the next-to-leading order perturbative 
QCD 0II3- The decay moments Z hv are calculated 
by using the charmed-hadron decay distributions given 
in Refs. [H[l6|. 



FIG. 3: Comparison of the total galactic plane tau neutrino 
flux with the total downward going atmospheric neutrino flux 
for the best fit values of the neutrino mixing parameters as a 
function of the neutrino energy. The intrinsic and oscillated 
parts of the total downward going atmospheric tau neutrino 
flux are also shown as a function of the neutrino energy. 



We show the relative contributions by the 7r, the K, 
and the charm decays to the overall flux in Fig. [21 
This is an extension of the Fig. 3 in Ref. [8j, where 
only the 7r and the K contributions are compared. It 
is obvious that the 7r decay contribution dominates for 
1 GeV < E < 10 GeV, while the K decay contribution 
dominates between 10 3 GeV and 10 5 GeV. The fraction 
of the charm-decay contribution rises rapidly at E > 10 5 
GeV and becomes dominant for E > 10 6 GeV. In this 
energy range, both the 7r and the K lose large fractions 
of their energies before decaying into the neutrinos. 

Additionally, the intrinsic atmospheric v T flux is also 
required to completely determine the total atmospheric 
v r flux. This flux is calculated using the perturbative 
QCD, since r neutrino arises from the D s decays jl7j . 
The flux is written as Eq. Q with Z p h replaced by the 
ZpD s and Zh VT replaced by the Zd sVt . We note that the 
Zd s v t contains two contributions. One arises from the 
decay D s —> v t t; the other follows from the subsequent 
tau-lepton decay, r — > v T + Y . The latter contribution 
is calculated using the decay distributions of the decay 
modes r — > v -rP, t —> v t tt, t — > v T a\ [rH lisj . and the 
t — > v T lvi |jj 0] (Fig. n shows the intrinsic downward 
going atmospheric muon and tau neutrino fluxes for 1 
GeV < E < 10 3 GeV). 

We have calculated the (^(E) by applying Eq. @ 
with (j) Vfi ^(E) given by d 2 N Vii T (E, X)/dEdX and inte- 
grating over the slant depth X. For £ < 70°, the os- 



cillation probability P(v^ 



is calculated using the 
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FIG. 4: The comparison of the galactic-plane and the down- 
ward going atmospheric v T fluxes in the presence of neutrino 
oscillations with maximal mixing as a function of the neutrino 
energy in GeV. For Sm 2 — 2.4- 10 -3 eV 2 , both fluxes cross at 
E = 2.3 GeV. 
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FIG. 5: The comparison of the galactic-plane v T flux and the 
atmospheric v T flux coming in the zenith angle £ = 180° in 
the presence of neutrino oscillations with maximal mixing as 
a function of the neutrino energy in GeV. For 8m 2 = 2.4- 10 -3 
eV 2 , both fluxes cross at E = 800 GeV. 



relation X = Xq exp(— L cos / ho) / cos£ with Xq = 1030 
g/cm 2 and ho — 6.4 km. Here L is the linear distance 
from the neutrino production point to the detector on 
the Earth 0. This gives, for instance, P{v^ — > v T ) ~ 
sin 2 (4.5 • 10" 2 (GeV/£0) at £ = 0°, for the best-fit values 
of the neutrino mixing parameters. 



IV. THE COMPARISON 

The comparison of the galactic-plane and the down- 
ward going atmospheric v T flux is given in Fig. and 
Fig. 01 in the presence of neutrino oscillations (we plot 
diV„/d(log 10 J5) instead of &N V /&E). The former flux 
clearly dominates the latter for E > 10 GeV, whereas 
the two fluxes cross at E = 2.3 GeV for Sm 2 = 2.4 ■ 10~ 3 
eV 2 and sin 2 2$ = 1. This comparison is however subject 
to the uncertainty of galactic-plane v T flux by the choices 
of the density n p and the distance R mentioned before. 

One can see that the atmospheric v T flux is sensitive 
to the value of 8m 2 for E < 20 GeV. Furthermore, a 
change of slope occurs for the atmospheric v T flux at 
E w 20 GeV. Beyond this energy, the slope of the atmo- 
spheric v T flux is identical to that of the galactic-plane v T 
flux. For E > 20 GeV, the atmospheric v r flux is 
intrinsic, i.e., coming from the D s decays, whereas the 
galactic-plane v T flux arises from the oscillation of the v^, 
which is produced mainly by the n decays. In both cases, 
the hadrons decay before interacting with the medium. 
Such a feature dictates the slope of the outgoing neu- 



trino flux. Below 20 GeV, however, the atmospheric v T 
flux predominantly comes from the oscillations, i.e., 
$£{E) re <j) Vii {E)-sm 2 29- sin 2 (L/L osc ) following Eq. ©■ 
Since L osc = 4E/6m 2 pa 330 km for E = 1 GeV with 
Sm 2 = 2.4- 10~ 3 eV 2 , we approximate sin 2 (L/L osc ) with 
(L/L osc ) 2 so that 4%?(E) - cj)^{E)E- 2 . Because the 
neutrino oscillation effect steepens the 0„ T spectrum for 
E < 20 GeV, the slope change of 0„ T at E ps 20 GeV is 
significant. 

We have also worked out the comparison of the 
galactic-plane and the atmospheric v T fluxes for sev- 
eral other zenith angles. For instance, the crossing of 
the galactic-plane and the atmospheric v T fluxes occurs 
at E = 6.0 GeV for the zenith angle £ = 60° for 
Sm 2 = 2.4 • 10 -3 eV 2 with the maximal mixing. Essen- 
tially, the atmospheric v T flux increases with the zenith 
angle £ while the galactic-plane v T flux stays unchanged. 
Moreover, from the phenomenological point of view, a 
larger £ only implies a larger atmospheric tau neutrino 
background while our focus is on the prospective obser- 
vation of the astrophysical tau neutrinos. 

Two key factors determine the angular behavior of the 
former flux. First, the atmosphere depth is larger for £ = 
60° compared to the downward direction. Second, the 
atmospheric muon neutrinos are produced more far away 
from the ground detector in this angle, implying a larger 
— > v T oscillation probability. We have found that, 
for E = 10 GeV, the downward going muon neutrinos 
are produced on average about 13 km from the ground 
detector. At £ = 60°, the above distance increases to 34 
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FIG. 6: An illustrative comparison of the downward going 
atmospheric and v T fluxes and the corresponding galactic- 
plane neutrino fluxes in the presence of neutrino oscillations as 
a function of the neutrino energy in GeV. The best fit values 
of the neutrino mixing parameters, namely, 5m 2 = 2.4 ■ 10~ 3 
eV 2 and sin 29 = 1, are used. The galactic-plane and the 
atmospheric fluxes cross at E = 5 • 10 GeV. 

km. In the angular range £ < 70°, the curvature of the 
Earth can be neglected and the angular dependence of 
the neutrino flux can be easily calculated |l9| . 

In Fig. [SI we show the comparison for the zenith an- 
gle £ = 180° (vertically upward going direction). For 
5m 2 = 2.4- 10~ 3 eV 2 with the maximal mixing the cross- 
ing between the two fluxes cross at 800 GeV. The oscilla- 
tory behavior in the total atmospheric tau neutrino flux 
is a manifestation of the narrow zenith angle dependence. 

We have so far neglected the 3-body muon decay con- 
tribution to the atmospheric flux. As mentioned ear- 
lier, this contribution is relevant only for 1 GeV < E < 
10 GeV. The ratio of this contribution to the total atmo- 
spheric flux is available in the literature 0, . For 
cos £ = 0.4 (£ = 66°), E = 1 GeV, the ratio is about 50%. 
It drops to 30% at E = 10 GeV for the same £ 16]. This 
ratio is not very sensitive to the zenith angle £. The er- 
ror due to neglecting this contribution propagates to the 
determination of the atmospheric v T flux by the neutrino 
oscillation effect. Numerically, this error is comparable 
to the uncertainty of the v T flux due to the uncertainty 
of the Sm 2 , as (E) w ^ (E) ■ sin 2 29 ■ {8m 2 ■ L/AE) 2 . 

V. CONCLUSIONS 

The results presented in Figs. 0] [S] and [S] indicate the 
opportunities for the tau neutrino astronomy in the GeV 



energies for the incident zenith angles 0° < £ < 180° 
in the two neutrino flavor mixing approximation. The 
galactic-plane tau neutrino flux dominates over the at- 
mospheric tau neutrino background until a few GeV's. 
Furthermore, for E < 20 GeV, the former flux has a 
rather different slope from that of the latter. This is 
an important criterion for distinguishing the two fluxes, 
particularly given that the normalization of the galactic- 
plane tau neutrino flux is still uncertain. 



We have pointed out that the dominance of the 
galactic-plane tau neutrino flux over its atmospheric 
background in GeV energies is unique among all the 
considered neutrino flavors. This is depicted in Fig. 
Because of the v ^ — > v T neutrino oscillations, the total 
galactic v T flux is identical to that of the galactic flux. 
However, the atmospheric flux is much larger than 
the atmospheric v r flux. As a result, in the presence 
of neutrino oscillations, the crossing energy value for the 
galactic-plane and the atmospheric fluxes is pushed 
up to 5 • 10 5 GeV, which is drastically different from the 
tau neutrino case. 



We have estimated the galactic-plane tau neutrino in- 
duced shower event rate N Vt for the forthcoming one 
Mega ton class of detectors. It is obtained by convolving 
the total galactic-plane tau neutrino flux shown in Fig. 
with the charged-current v T N interaction cross section 
(using the CTEQ6 parton distribution functions [2(]])- It 
is found that the N Vt can be 4.5 (7.5) for E > 10 GeV 
with a data taking period of 3 (5) years. 



A remark concerning the possible background induced 
by the electron neutrino events and the neutral current 
events to the prospective observation of tau neutrino 
events is in order. Below 10 3 GeV, the tau lepton decay 
length is less than a mm. This tau lepton is produced in 
the detector in the galactic-plane tau neutrino induced 
interactions. There are certain specific signatures of the 
tau neutrino induced tau leptons such as the appearance 
of the kink at the tau lepton decay (absent for the elec- 
trons) 21], as well as the relative characteristic fractional 
energy sharing from the incident neutrino |22j| . A large 
scale finely grained detector with a resolution of a few /xm 
will thus be required to disentangle the galactic-plane tau 
neutrino induced events from the events induced by the 
electron neutrinos and the neutral current events on the 
event-by-event basis. 
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